source LoadPDB.tcl 

buildGroupList step5_assembly.pdb {N C1? O1* P } by residue 

puts [llength $groups]
puts $groups 
